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Abstract 

Recent Chandra X-ray observations of jets in central cluster galaxies show interesting 
features around the symmetry plane. We have carried out simulations involving 
bipolar jets, removing the artificial boundary condition at the symmetry plane. We 
use a very low jet density (IGM/jet ~ 10 4 ) and take into account a decreasing 
density profile. We find that the jet bow shock undergoes two phases: First a nearly 
spherical one and second the well-known cigar-shaped one. We propose Cygnus A 
to be in a transition phase, clear signs from both phases. Due to inward growing 
of Kelvin Helmholtz instabilities (KHI) between cocoon and shocked IGM, mass 
entrainment is observed predominantly in the symmetry plane. We propose this 
mechanism to produce some of the so far enigmatic X-ray features in the symmetry 
plane in Cygnus A. 
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1 Introduction 



Radio galaxies are traditionally 
known by their radio emission, which 
arises from magnetised jet plasma, 
accelerated to large velocities in the 
active nucleus, and expelled up to 
Mpc distances. During their growth, 
they displace and compress the ambi- 
ent gas. In the central and dominant 
cluster galaxies (cDGs), this gas is 
usually dense enough to be seen in 
the X-ray regime by bremsstrahlung. 
The Chandra observatory can now 
resolve these centers of galaxy clus- 
ters with unprecedented resolution. 



The displaced hot gas (10 7 K) has 
been detected in several cases (e.g. 



ISoker et"aL 120021 and references 
therein). An exceptional case is the 
Cygnus cluster with its famous FR II 
radio galaxy Cygnus A. This is the 
only nearby (redshift z=0.06) pow- 
erful jet source. In Cygnus A, the 
leading bow shoc k can be found on 



the X-ray image (|Smith et all l2002i 



also compare Chandra homepage). 
The identification is easy: It closely 
engulfs and follows the radio lobes in 
vicinity of the hotspots, and forms 
an almost spherical bubble around 
the radio emission. It also reproduces 
the armlength asymmetry. From this 
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data, the axis ratio of the bow shock 
can be determined to 1.2. Compar- 
ison with the radio data gives a 
bow-shock-to-cocoon width of 3.8. 
In the symmetry plane, X-ray bright 
filaments are present, next to and 
within the radio lobe. Such excellent 
data demands detailed modelling. 
We have simulated a jet with the 
conditions assumed to be present in 
Cygnus A (sec. 2). The simulation is 
bipolar in order to take into account 
interactions of the backflows of the 
jets. Because it fails to reproduce 
the data, quantitatively, we present 
a parameter study in sec. 3. Finally, 
we present a modified jet model for 
Cygnus A in sec. 4 and discuss the ap- 
plication to comparable jet sources. 



2 Bipolar 3D Simulation 




Fig. 1. Density slice for the 3D simula- 
tion. The jet is injected in the center in 
opposite directions. 

extensions. The aspect ratio of the 
bow shock is monotonically increas- 
ing for the whole simulation. At the 
time shown, one jet arm is longer 
than the other one by « 10%. The 
background is slightly asymmetri- 
cally perturbed, with the higher den- 
sity on the side of the longer jet. 
This shows that the asymmetry is 
caused by the action of KHIs on the 
jet beam. The central part of the 
bow shock is pressure driven and its 
width turns out to scale as t 0,6 , for 
the whole simulation time. 



Using the code Nirvana (IZiegl er and York 
19971 ?. we simulated dKrausel . I2002T T" 
a bipolar jet in 3D and cylindri- 
cal coordinates with the param- 
eters of the Cygnus A je t , est i- 
mated by Carilli and Barthel ( 1996j ) : 
beam radius rj = 550 pc, velocity 
Uj = 0.4c, power Lj = 10 46 erg/s, 
and external density profile p(r) = 
10~ 25 [1 + (r/35kpc) 2 ]~ 9 / 8 g/cm 3 . 
The density contrast t] is 6 x 10~ 3 
at the inlet and the Mach number is 
M = 10. 

A slice through the computational 
volume at a time of 1.6 Myr is shown 
in Fig. 1. Some basic features of such 
a jet are: Initially, the bow shock 
forms a spherical bubble. Then it 
first elongates the bubble, and after 
some time it produces cigar shaped 



-The observed jet in Cygnus A has ba- 
sically the same features as the sim- 
ulated one. The bow shock is a big 
bubble with a cigar shaped extension 
on one side. The armlength asymme- 
try is ~ 10%. But the bow-shock-to- 
cocoon width £z 1.3 and the axis ratio 
of the bow shock (> 3) is far-off the 
observed values. Hence, the assumed 
parameters wer not a good choice. 



3 Parameter Study 



In order to find out, what parame- 
ters differ and in which direction, we 
performed an axisymmetric study 
with constant background, density 
contrasts rj G [10~ 5 ; 10~ 2 ], a rid Mach 
numbers of 3, 8, and 26 ( Krausel . 
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stant density profile, one obtain s 
the solution of ICastor et al. ( 1975| ): 
r = (5L/47rp )°' 2 t a6 - The criti- 
cal radius for the first deviation 
from spherical symmetry is given 
by: fi/rj = 0.5?7 _1 / 4 which can be 
checked in the simulations. All the 
simulated jets show a high cocoon 
pressure, which decreases in time, as 
expected for a blastwave with con- 
stant energy input. The jet beam 
does not react by compression or 
expansion, but by adjusting its own 
pressure via oblique shock waves. 
The result is approximate pressure 
equilibrium, besides the hot spot. 
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Fig. 2. Logarithmic density plots of the 
parameter study. Top row: rj = 10 -2 , 
bottom row: T] = 10~ 4 . Left column: 
M=3, right column: M=26. Darker 
colour corresponds to less dense regions. 



2003). Fig. 2 shows that the trend 
foun d in classical p arameter stud- 
ies f Normanl . ll993L cocoon radius 



r c oc M77 -1 / 4 ) is also confirmed here: 
the lighter jets show fatter cocoons. 
However, the cocoons of the 77 = 10~ 4 
jets are still under-expanded, which 
may be the reason for the indepen- 
dence of r c on M. The bow shocks are 
spherical at the beginning, because 
the pressure inside the jet bubble 
drives them faster than the thrust at 
the jet head. The relevant equation 
is the force balance at the bow shock, 
which can be transformed into the 
following integral equation 

r t ti 

J M(r)rdr = 2 / dt x J E(t 2 )dt 2 (1) 



that relates the mass M inside 
the radius r to the energy E, in- 
jected during the time t. For con- 
stant energy injection rate and con- 



4 Discussion 



The 3D simulation has falsified the 
assumed parameters for Cygnus A. 
The observed cocoon and bow shock 
is wider than in the simulation. This 
points to a more extreme density con- 
trast. It turns out from the param- 
eter study that 77 ps 10~ 4 is needed 
in order to get the observed values. 
This ag rees well with the e xtrapola- 
tions of Rosen et al. ( 19991 ) concern- 
ing the cocoon width. Equation 1 
can be used to infer the average jet 
power (L = 8 x 10 46 erg/s), and the 
age (t = 27 Myr), which agrees with 
estimates in the literature. The simu- 
lations have shown that the pressure 
has the same order of magnitude in 
the whole system. The external pres- 



sure i n the cluster gas is ((Smith et al 
20021 ): p cxt ^3x 10~ 10 erg/cm 3 . The 



magnetic energy density in the hot 
spots is Wm ^us = 9 x 10 ~ 10 erg /cm 3 



or higher (j Wilson et all |2 



as- 
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suming self-synchrotron-Compton 
origin of the X-ray emission. This 
is roughly the equipartition value. 
Since the magnetic field in the hot 
spots is probably shock- amplified, 
this implies a similar value of the 
internal energy in the jet and the 
cluster g predicted by the 

above model. Given the density con- 
trast of rj = 10 -4 and taking into 
account a possibly relativi stic jet 
( compare IScheck et al. . 2002j ). we es- 



timate: r 2 /ipjc 2 « 10~ 8 erg/cm 3 (r: 
bulk Lorentz factor, h: specific rela- 
tivistic enthalpy, pf mass density). 
Consequently, the kinetic power of 
the jet is: L^ in = 2nr?pjT(Th — 
l)pc 3 = 5 x 10 45 (1 - 1/Th) erg/s. 
(P = Vj/c) Since this falls short 
of the value derived above, the en- 
ergy has to be in the magnetic 
power: L B = r 2 T 2 f3cB 2 / '4 = 5 x 



10 



1 1 , 



magj/ ^mag,HS J 



\2p2 



It follows, 



that T has to be roughly 20. The 
mass flux in the jet is then given by: 
M = irrfp;T0c = 3x l Q -^A Wyr. 
Appl and Camenzind ( 1988| ) esti- 



mate the mass flux needed in or- 
der to produce the observed lumi- 
nosity in the Cygnus A hot spots 
downstream of a MHD shock: M = 
5 x 10~ 4 e~jM Q /yr (e re i: energy frac- 
tion in relativistic electrons). This 
does not preclude protons in the jet, 
and places h in the regime of a few. 

The 3D simulation shows that the 
KHI is strongest in the symmetry 
plane, drawing fingers of shocked 
external cluster gas into the radio 
cocoon. We propose that some of 
the X-ray enhancements inside the 
cocoon could be due to this process. 
Outer X-ray enhancements may be 
caused by waves that the KHIs excite 



in the shocked cluster gas. 

Some cDGs show similar features like 
Cygnus A. For exam ple, the X-ray 



structure of 3C 317 (IBlanton et al 



(B 



20mh (two elliptical rings, elongated 
in the same direction) can be ex- 
plained bubble with a cigar 
shaped extension, seen at an inclina- 
tion of 37°. 
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